Excitonic spin density wave state in iron pnictides 
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We examine the appearance of a spin density wave in the FeAs parent compounds due to an 
excitonic instabihty. Using a reahstic four-band model, we show that the magnetic state depends 
very sensitively upon the details of the band structure. We demonstrate that an orthorhombic 
distortion of the crystal enhances the stability of the antiferromagnetic order. 
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Introduction. The superconductivity of materials con- 
taining FeAs layers is currently receiving much attention. 
Like the cuprates, these systems become superconducting 
upon chemical doping of an antiferromagnetic (AF) par- 
ent compound, specifically ReFeAsO (Rc is a rare earth 
ion) or AeFe2As2 (Ae is an alkaline earth ion).^-^ Intrigu- 
ingly, the AF state occurs only in the presence of an 
orthorhombic distortion of the crystal, which fixes the 
AF ordering directioni^ The likely role of spin fluctu- 
ations in producing the superconductivity has lead to 
intense scrutiny of the AF phase. The relatively small 
value of the moment at Fe sites (^i^i^ metallic transport 
properties f^i^ and observations of reconstructed Fermi 
surfaces,—!^ provide strong evidence that the AF state 
is a spin density wave (SDW) arising from the nesting of 
electron and hole Fermi surfaces 

In analogy to Cr^i^ a theory of the SDW based 
upon the excitonic pairing of electrons and holes has 
been proposed.— i^^ii^ It is important to determine if 
this scenario is sufficient to explain the AF state, or 
whether a more complicated multi-orbital approach is re- 
quiredi^'i^ii^'^ii^ As previous works have used a highly- 
idealized model of the electronic structurC j^'^i^'*'^^ with 
only two Fermi surfaces instead of the likely four or 
more)2ii£ it is not clear if the excitonic SDW can give the 
observed magnetic ordering. '^'^ Furthermore, the effect of 
the orthorhombic distortion on such a state remains un- 
known. We address these problems here by studying the 
appearance of the excitonic SDW in a four band model 
of LaFeAsO.^^ Using a mean- field theory, we show that 
the SDW state is sensitively dependent upon the doping 
and the details of the band structure.-- In particular, we 
examine the response of the SDW phase to changes in 
the ellipticity of the electron pockets, the relative size of 
the hole pockets, and an orthorhombic distortion of the 
crystal. 

Theoretical model. We model the FeAs planes as a 
2D interacting four-band system where two bands have 
electron-like Fermi surfaces and the other two have hole- 
like Fermi surfaces. We write the Hamiltonian as 

n=l,2 k,cr 
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where cj^^.^ (/Ike) creates a spin-cr electron with momen- 
tum k in the electron-like (hole-like) band n. Due to 
the out-of-plane arrangement of the As ions, the crys- 
tallographic unit cell of the FeAs plane contains two 
Fe ions. Our band structure is given in terms of this 
unit cell, but in the discussion of magnetic properties 
it is more useful to refer only to the Fe lattice, which 
requires us to "unfold" the Brillouin zonci^ Assuming 
crystallographic unit-cell dimensions a x a, the bands 
with electron-like Fermi surface have dispersion ej^j^ — 
ee + te,i[cos{tj,a) + cos{kya)]+te,2 cos{[kx + {-iy''ky]a/2), 
while for the hole-like bands we have e^^. — eii,n + 
t/i,ria[cos(fca;a) -I- cos{kya)] + tu.na cos{kxa) cos{kya). In 
units of eV, we use — 1.544, t^.i = 1.0, ie,2 = —0.2, 
eh,i = -0.335, thxi = 0.24, th,ia = 0.03, = -0.512, 
th,2,i — 0.315, and th,2,2 = 0.06. We keep only the 
bands which intersect the Fermi surface. For electron 
filling Tici = 4, corresponding to the undoped parent 
compounds, we find the dispersion and Fermi surface as 
shown in Fig.s[IJa) and (b), respectively. Note that the 
nesting of the hole and electron Fermi surfaces is not per- 
fect, since both the shape and the enclosed area differ. 
Our model reproduces the Fermi surface and low-energy 
velocities of the band structure proposed in Ref. [ll| for 
LaFeAsO, but unlike Ref. I 111 obeys the correct periodicity 
of the Brillouin zone. 

The interaction terms in Eq. ([T]) describe a density- 
density interaction and correlated transitions between 
the electron and hole bands, with contact potentials 
gi > and g2 > 0, respectively. At low temper- 
atures, the system is unstable against the pairing of 
electrons and holes, producing an excitonic stateJ^i^iiS^ 
Although a rich variety of excitonic phases are possi- 
ble, here the SDW state has the largest effective cou- 
pling constant gs = gi + 2^2'^'^ At mean- field level, 
we therefore decouple the interaction terms via the in- 
troduction of the real SDW excitonic averages Aeh = 
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FIG. 1: (color online) (a) Dispersion of the electron and hole 
bands along high-symmetry directions of the Brillouin zone, 
(b) Fermi surface at filling n^i — 4. The bands are represented 
by the same curves as in (a), (c) Variation of the excitonic 
gaps with temperature at n^i — 4. 



FIG. 2: (color online) (a) Electron Fermi surfaces at St = 0.2 
(thick lines) compared to 5t — (thin lines), (b) Variation 
of Ae,h with 5t at T = IK. (c) Variation of Ae,i with St at 
various temperatures. Shading beneath the curve indicates a 
single-Q solution show in (b) and (c). 



i9s/V)T,kJ2a <^i4M+Q,afhka) where Q = (7r/a,7r/a) is 
the nesting vector [see Fig-IHb)] and e{h) takes values of 
1 or 2 to index the electron (hole) bands. A^h is regarded 
as the order parameter of the SDW stateji^i^i although 
it is only indirectly related to the staggered magnetiza- 
tioni^^ As each electron pocket is mapped to a different X 
point of the enlarged Brillouin zone upon unfolding,^ the 
Aih and involve orthogonal nesting vectors Qi and 
Q2 with respect to the Fe sites, respectively. When both 
Aih and A2/1 are non-zero, therefore, the magnetization 
is the superposition of two orthogonal SDW states, each 
with stripe-like orderingJi 

After decoupling the interaction terms, we obtain the 
equilibrium mean-field solution by numerical minimiza- 
tion of the free energy F with respect to the A^.h- This 
was calculated over the 2D Brillouin zone with at least a 
1000 X 1000 k-point mesh. Throughout this work we set 
the effective SDW couphng constant to be = 0.84eV, 
as at Uei = 4 this gives a partially-gapped Fermi surface 
in the SDW state with reasonable critical temperatures: 
as shown in Fig. [Ijc) we find that A^i is non-zero below 
TsDWi — lOOK, while Ae2 appears below Tsdw2 — 6.5K. 
When all four averages are non-zero, we find the inequal- 
ity A11A12A21A22 < 0; when only two A^h are present, 
their signs are independent. 

Ellipticity of the electron pockets. As seen in Fig.[ljc), 
both electron bands participate in the excitonic insta- 
bility at rici = 4. This corresponds to a Qi + Q2 
SDW, whereas only a single-Q SDW is experimentally 
observed.^ It has previously been noted that these two 
SDW phases should lie at similar energiesr^ and so it 
is interesting to see whether slight changes in the band 
structure can stabilize a single-Q state. This might 



be achieved, for example, by reducing the ellipticity of 
the two electron pockets so as to enhance their com- 
petition for the same states in each hole band. We 
therefore modify the electron dispersions e^j,. ej^j]^. -I- 
2(— l)"ie,2i5i sin(fca;a/2) sin(A;j,a/2), where the dimension- 
less parameter St controls the ellipticity of the electron 
pockets. We compare the electron pockets at St — 0.2 
and (Si = in Fig.[2i;a). 

We find that even very small values oi St can qual- 
itatively alter the mean-field state. The evolution of the 
Aeh with increasing St at T — IK is plotted in Fig.[2I^b). 
Reducing the ellipticity of the electron pockets tends to 
suppress the excitonic state, with Ae2 disappearing be- 
fore St = 0.01 is reached. At St = 0.108 the system 
undergoes a first-order transition from the Qi + Q2 state 
into a single-Q state. A single-Q state is hence possi- 
ble at mean-field level by subtle modification of the band 
structure. Note that the single-Q states with nesting 
vector Qi and Q2 are degenerate. Further increasing 
St, the system undergoes a first-order transition into the 
nonmagnetic state at St « 0.145. 

The variation of Ae.i with St at higher temperature is 
shown in Fig. [U^c); in all cases Ae.2 = 0. The first-order 
transition from the Qi + Q2 into the single-Q state only 
survives up to T 30K; at higher temperatures, the 
nonmagnetic state is reached from the Qi -I- Q2 phase by 
a second-order transition. Interestingly, we see that the 
critical value of St increases with T, even as the value of 
Aei at St = is suppressed. This re-entrant behaviour 
is a generic feature of the phase diagram of the excitonic 
insulator , ^^'^^ and may indicate the presence of a low-T 
incommensurate SDW stateJ^^ 

Hole pocket disparity. The SDW state is sensitively 
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FIG. 3: (color online) (a) Hole Fermi surfaces at 5h = 0.05eV 
(thick lines) compared to Sh — (thin lines) at rid = 4. (b) 
Dependence of the critical temperature Tsdw of the excitonic 
phases on electron filling n^i for different values of Sh. (c) 
Variation of the Aeh with rioi at T = IK and 5h — 0. (d) 
Same as (c) but for Sh = O.OSeV. 



dependent not only upon the shape, but also upon the 
size of the Fermi surfaces. This can be demonstrated 
in two ways: by raising the energy of the second hole 
band ^ — * ^2 k + ^'^ 'Otmt the two hole Fermi sur- 
faces converge together, or by varying the fiUing Uoi to 
improve the nesting between one of the hole Fermi sur- 
faces and the two electron pockets. At bh = O.OSeV, 
the two hole Fermi surfaces are nearly coincident when 
Uei = 4, see Fig.lSJa). As shown in Fig. [3Jb), this energy 
shift strongly alters the nd-dependence of the maximum 
temperature Tsdv^ at which at least one Aeh non-zero. 
When 5h, = 0, our model displays two distinct peaks 
in the Tsdw vs. riei curve, with a sharp minimum at 
Uoi ~ 3.99. This behaviour qualitatively disagrees with 
experiment, which shows only monotonic suppression of 
Tsdw with electron-doping.^ The behaviour of TgDW at 
bh — 0.05eV is in much better agreement with experi- 
ment, with only a single maximum. Note that the maxi- 
mum value of Tsdw in both cases is comparable to that 
in the ReFeAsO systems<^ 

The Tsdw vs. n^x curves can be understood by ex- 
amining the evolution of the Ag/j with nd at T = IK, 
plotted in Fig. [HJc) for = and in Fig. ^6) for 
(5/i = 0.05eV. Note that the values of the pairs (An, A22) 
and (A2i,Ai2) may be swapped at every point. At 
6h — 0, the two distinct peaks in Fig. ^}a) correspond to 
a maximum in |Aei| for electron doping and in |Ae2| for 



hole doping. The maximum values are different due to 
different densities of states in the hole bands. These max- 
ima occur when the area enclosed by the hole Fermi sur- 
face is the same as that enclosed by each electron Fermi 
surface. It is interesting to note that at both maxima a 
single-Q state is stable. 

When bh — O.OSeV, the conditions for |Aei| and |Ae2| 
to display a maxima coincide at rioi = 4, as the area en- 
closed by each hole Fermi surface is almost equal. We 
hence see a complicated coexistence between the four or- 
der parameters: at weak hole doping, all four Ag^^ are 
non-zero; at weak electron doping, the excitonic instabil- 
ity of the two electron bands involve different hole bands. 
Although Aei is dominant over most of the doping range, 
at extreme hole doping a state with only Ae2 non-zero is 
realized, corresponding to the weak asymmetry seen in 
the Tsdw vs. nd curve in Fig.[3Ib). 

Orthorhombic distortion. In all known FeAs parent 
compounds, the SDW phase occurs only in the presence 
of an orthorhombic distortion of the crystal. It is found 
that the stripe- like SDW has its nesting vector Q oriented 
along the longer crystal axisi^i^ Here we see how this can 
be understood within our model on the basis of the effect 
of the orthorhombic distortion on the Fermi surfaces. 

Under an orthorhombic distortion, the energy shift of 
a state with wave- vector K in the unfolded Brillouin zone 
is Se-K ~ J2a p KaYa,f3Kp where Ya^p is the strain ten- 
sor and we have Y^x — —Yyy and Y^y = 0^ Note that 
the wave-vectors K in the unfolded Brillouin zone are 
rotated by 4S° with respect to the wave-vectors in the 
crystallographic Brillouin zone. We approximate the en- 
ergy shifts Se-K by their value near the chemical poten- 
tial, as the Fermi surface shape dominates the physics of 
our model. The energy shifts of the hole states near the 
zone centre are therefore neglected, as they will be much 
smaller than those experienced by the electron pockets. 
Since the electron pockets are small, we assume that their 
energy shifts are isotropic. Furthermore, the sign of the 
energy shift will be opposite for the electron pockets at 
the X points along the axes of compression (negative 
energy shift) and dilation (positive energy shift).— We 
hence model the effect of the orthorhombic distortion by 
ej^^ — > e^j^ + (— l)"(5e. We compare the electron pockets 
at Se = 0.04eV and ,5e = in Fig. g^a). 

The dependence of the A^h upon Se at T ~ IK and 
T = 40K is plotted in Fig. gl^b) and Fig. |D^c), respec- 
tively. The effect oi de ^ is to enhance the pairing be- 
tween the larger electron and hole pockets (An) and also 
the smaller electron and hole pockets (A22), while sup- 
pressing the pairing between the smaller electron (hole) 
and larger hole (electron) pockets. In analogy to the ef- 
fect of doping, this can be readily understood as due to 
the changes in the area enclosed by each electron Fermi 
surface. Due to the enhanced excitonic pairing, the free 
energy F shows monotonic decrease with increasing Se, 
see Fig. IHJd) . As the orthorhombic distortion should in- 
crease the elastic energy of the lattice, it is therefore pos- 
sible that the total free energy of the crystal will show a 
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FIG. 4: (color online) (a) Electron Fermi surface pockets at 
Se = 0.04eV (thick lines) compared to (5e = (thin lines), (b) 
Dependence of the A^h upon 5e a.t T = IK. (c) Same as for 
(b) but at T = 40K. (d) Difference AF = F - Fo per site 
between the free energy F and its value in the normal state 
Fo at Se = 0. 



minimum at a non-zero value of the distortion. Deeper 
investigation of this scenario is left for future work. 

The T = 40K case shows a large range of Se where An 
is the only non-zero excitonic average, i.e. the distor- 
tion stabilizes a single-Q SDW state due to the enhanced 
nesting between the larger electron and hole pockets. In 
contradiction to experiment, however, the Q vector is ori- 
ented along the shorter crystal axis. This does not neces- 



sarily invalidate the excitonic scenario: our model Eq. ([T]) 
has equal coupling constants between the different bands. 
Were hole band 2 to interact more strongly with the elec- 
tron bands than hole band 1, so that |Ae2| ^ l^eil in 
the undistorted system, the enhancement (suppression) 
of A22 (A12) by the orthorhombic distortion would likely 
stabilize a SDW state with the observed Q vector. 

Conclusions. We have presented a mean-field study of 
the excitonic SDW state for a realistic four-band model 
of the FeAs parent compounds. We find that the SDW 
state is sensitively dependent upon the band structure. 
For a tetragonal unit cell, a two-Q SDW is realized at 
rici = 4; small changes in the electron pocket ellipticity 
or the doping, however, stabilize the observed single-Q 
state. Varying the relative size of the hole pockets qual- 
itatively changes the Tsdw vs. rid curve, agreeing best 
with experiment when the hole pockets are almost coin- 
cident.^''^ The dominant effect of an orthorhombic distor- 
tion of the crystal on the band structure was identified as 
altering the size of the electron pockets. This changes the 
nesting condition between the Fermi surfaces, and can 
realize a single-Q SDW. Our analysis suggests that the 
electron pockets interact more strongly with the smaller 
hole Fermi surface than with the larger. We conclude 
that the excitonic SDW model is capable of qualitatively 
describing the AF phase of the FeAs parent compounds. 
The strong sensitivity of the SDW state upon the band 
structure, however, shows that a quantitative description 
requires a more detailed understanding of the electronic 
structure than is currently available. 
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